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Chain inflation proceeds through a series of first order phase transitions, which can release con- 
siderable gravitational waves (GW). We demonstrate that bubble collisions can leave an observable 
signature for future high-frequency probes of GWs, such as advanced LIGO, LISA and BBO. A 
"smoking gun" for chain inflation would be wiggles in the spectrum (and consequently in the tensor 
spectral index) due to the multiple phase transitions. The spectrum could also be distinguished 
from a single first order phase transition by a small difference in the amplitude at low frequency. 
A second origin of GWs in chain inflation are tensor modes from quantum fluctuations; these GW 
can dominate and be observed on large scales. The consistency relation between scalar and tensor 
modes is different for chain inflation than for standard rolling models and is testable by Cosmic 
Microwave Background experiments. If inflation happened through a series of rapid tunnelings in 
the string landscape, future high frequency probes of GW can shed light on the structure of the 
landscape. 
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Several experiments are underway to detect the 
stochastic gravitational background from early universe 
cosmology. These gravity waves (GWs) propagate almost 
freely throughout the entire history of the universe and 
thus can be a direct source of information about the uni- 
verse at very early times. GW experiments include the 
ground-based LIGO [l| and VIRGO [2j taking data now, 
and the upcoming space-based LISA Q . In the more dis- 
tant future are the proposed space-interferometers, BBO 
[H, with possible launch within 20 or 30 years. 

One source of GWs is inflation @ , a quasi-exponential 
growth phase of the universe which can explain its homo- 
geneity and isotropy as well as the generation of density 
perturbations required for structure formation. There 
are two types of contributions to GW from inflation. 
In slowly rolling models, quantum fluctuations of space- 
time lead to contributions proportional to the height 
of the inflaton potential. In tunneling models, there is 
an additional contribution due to bubble collisions of 
true vacuum bubbles at the end of the phase transition 
, [ll^ ] . In this paper we consider the model of chain 
inflation, which gives rise to both types of contributions. 
Regarding inflation, on ~ horizon scales, WMAP con- 
strains the GW amplitude, tt G wh 2 < 10" 15 @. Future 
proposed cosmological probes, e.g., CMBPOL, would be 
able to detect GW from inflation at large scales, if the 
energy scale of inflation ~ GUT scale Q . 

In chain inflation Q, the universe tunnels rapidly 
through a series of first order phase transitions. One can 
imagine a multidimensional potential containing many 
minima of varying energy. The universe starts out in 
a high-energy minimum, and then sequentially tunnels 
down to bowls of ever lower energy until it reaches the 



'Electronic address: 
t Electronic address: 



amjad@umich.edu 



bottom. During the time spent in any one of these min- 
ima, the universe inflates by a fraction of an e-fold. After 
many hundreds of tunneling events, the universe has in- 
flated by the 60 (or so) e-folds required to resolve the 
cosmological problems. In one variant of chain inflation, 
the tunneling corresponds to quantized changes in four 
form fluxes [9( such as may be found in the string land- 
scape. Other variants have been discussed in [l(| Hlj . 
but see [l2j]. At each stage, the phase transition is rapid 
enough that bubbles of true vacuum intersect one another 
and percolation is complete. Thus the failure of "old" 
inflation [j| to reheat is avoided, and "graceful exit" is 
achieved. 

In the zero temperature limit, the nucleation rate T 
per unit spacetime volume at a phase transition for pro- 
ducing bubbles of true vacuum in the sea of false vacuum 
through quantum tunneling is T(t) — Ae~ SE , where Se 
is the Euclidean action for the bounce solution extrapo- 
lating between false and true vacua [l3[ and A is a deter- 
minantal factor. For a first order phase transition, with 
Einstein gravity, it has been shown that the probability 
of a point remaining in the false vacuum is given by 
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Writing equation |T| as exp(— t/r), the lifetime of the 
field in the false vacuum is r = j^jj- The number of 
e-foldings for the tunneling event is 
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7 has to be greater than critical value, 7 C , where [14[ 

7c - 9/4tt, (4) 
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to achieve percolation and thermalization. This corre- 
sponds to having an upper bound on the number of e- 
foldings that that can be obtained in each stage of infla- 
tion: 
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The number of e-foldings required for successful in- 
flation depends on the energy scale of inflation and the 
reheating temperature. For a model with energy scale ~ 
GUT scale, about 60 e-foldings are required. This cor- 
responds to having at least ~ 200 tunneling events for 
GUT scale chain inflation. 

A study of density perturbations from chain inflation 
has been done by [ToL [H[; the former set of authors 
found that the right perturbations to match data can be 
generated. In this letter we calculate the gravity spec- 
trum from chain inflation. Since chain inflation proceeds 
through first order phase transitions, one expects a con- 
siderable amount of GWs from chain inflation. 

The spectrum of gravity waves (GWs) from multi- 
bubble collisions at a single phase transition (PT) with 
energy difference e between false and true vacua was 
worked out numerically in 16]. In this paragraph we 
review their results, and in the next section generalize 
to the multiple PT of chain inflation. We note that an 
example of a single PT model which resolves the graceful 
exit problem of old inflation by having time-dependent 
nucleation rates is in [2(| , [Il[ . The authors of [l6[ found 
that the fraction of vacuum energy released into gravity 
waves has the following spectrum: 
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The peak of the spectrum is at 
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Redshifting the frequency and energy of gravitational ra- 
diation as a -1 and a~ 4 respectively, at the current epoch 
the peak frequency, /o, and the corresponding logarith- 
mic contribution of GWs to critical density at the peak 
are given by [l9| : 
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Here Qq is the GW energy density in units of the critical 
density p c required to close the universe and h is today's 
Hubble constant in units of 100 km/s/Mpc. T* is the 
temperature increase right after a single phase transition, 
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and the total number of relativistic degrees of freedom 
at temperature T* is taken to be ~ 100. We note 
that in the past few years, the GW resulting from bubble 
collisions have been reexamined by and [Hj|. The 
former find that the fall-off at the low-frequency tail ~ 
rather than f^ 1 ' 8 ; this result leads to a smaller fall- 
off of GW amplitude and hence improves detectability. 
To be conservative, we will mainly stick to the result of 
fl(| and compare with the consequences of [ijj as needed. 

Spectrum of GWs from Bubble Collisions: We 
now generalize these results to a series of PT to obtain 
GW from a series of bubble collisions in chain inflation. 
Let us assume that the total number of minima that the 
universe passes through during chain inflation is N. As 
we will see its precise value will not have any impact on 
our final result. Then we will have N — 1 PTs during the 
course of chain inflation. The energy difference between 
two consecutive minima is assumed to be the constant 
value of e. We also assume that 7 does not change and it 
satisfies the condition 7 > j c . Assuming instant reheat- 
ing after each PT, the temperature of the universe rises 
after each stage by an amount T*, where T* is given in 
equation (fT2")) . The universe's temperature redshifts by 
a factor of exp(— \) during each bout of chain inflation. 
Thus the temperature of the universe at the end of chain 
inflation is 

T e = T*(l + e-*+- • ■ + e- {N - 1 ^) = T« K n - > (13) 

(1-e 

For N\ S> 1, the final temperature of the universe is 
simplified to 
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As stated above, every first order PT can release a sub- 
stantial amount of GWs from bubble collisions whose 
profile, peak frequency and energy fraction at the peak 
frequency are given by eqs. (O, (J7]) and (HJ). After 
each PT energy density in the gravitational radiation de- 
creases like a -4 , whereas the frequency redshifts like a -1 
[26| The peak frequency and fraction of critical density 
from GWs at the peak frequency generated during the 
rn-th PT redshifts as: 
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FIG. 1: Solid line shows the gross features of the GW spec- 
trum from bubble collisions in chain inflation for e 1//4 = 
10 8 GeV and X = |i where e is energy difference between 
vacua and \ is number of e-folds per stage. For more detailed 
features near the peak, see Fig. 2. The dotted line shows the 
GW spectrum generated from a single PT whose peak fre- 
quency and amplitude at the peak coincide with those from 
chain inflation; the small difference in the amplitudes at the 
low frequency tail is detectable. 

where a m , a e and ao are respectively the scale factors at 
the end of m-th phase transition, at the end of chain in- 
flation and today. The scale factor ratio a m /a e is equal 
to exp [— (N — m — 1 ) x] ■ If we also assume that the evolu- 
tion from the end of chain inflation up to today has been 
adiabatic [i.e. S oc a 3 g(T)T 3 = const], we have 

a e 14 /100\ 1/3 /l GeV\ . , 

8x10 Ur) (— )■ (17 » 

Thus 

Qh 2 {f m ) = 1CT 6 (1 - e-*) y e -4(A- m -i )x 1/3 (J9) 

The spectrum at other frequencies get redshifted such 
that the profile given in eq.© is preserved: 

( fih 2 (/ m )(f- ) 2 ' 8 f<f m 

[n/i a (/ TO )(£) f>fm, 

The total energy density in GW due to bubble colli- 
sions in chain inflation is obtained by summing all the 
contributions from the multiple tunneling events, 

JV-l 

n huh h 2 (f) = ]T n m h 2 (f) (21) 

m— 1 



FIG. 2: Tensor spectral index vs. frequency (blow-up of Fig. 
1 near the peak frequency). For / <C / pca k iit = —2.8. The 
wiggly structure in the spectrum at / < / pca k is a unique 
signature for chain inflation. For / > / pca k, «t jumps to 
— 1.8. Alternatively, if we follow [17|, the tensor spectral index 
takes the value of —1 at the low- frequency tail. 

and is plotted in Fig. (TT]). The frequency range plotted 
is from 10 Hz (roughly the minimum frequency at- 
tainable by LISA and BBO) to 10 2 5 Hz (roughly the 
maximum sensitivity frequency of Advanced LIGO [22l]). 
In the figure, we take N — 1000 minima, e 1 / 4 = 10 8 
GeV, and x — Xc = 1/3 (the critical nucleation rate). 
For such parameters the tip of the spectrum is observ- 
able in Advanced LIGO which is sensitive within a fre- 
quency band of 1-1000 Hz to stochastic signals with 
Hqw'* 2 ^ 10~ n . Part of the low-frequency tail of the 
spectrum could be seen at BBO too, which could probe 
GWs with much smaller amplitude, i.e. Hq^H 2 > 10~ 17 
for / - (10~ 4 - 10 1 ) Hz. 

One can see that the GW spectrum due to bubble col- 
lisions at frequencies smaller than the peak frequency is 
blue. This is in contrast with slow-roll inflation which 
can not produce a blue gravity spectrum on an extended 
range of scales, unless one violates the null energy condi- 
tion H > 0. This is rooted in the sub- horizon mechanism 
of generating the GWs from bubble-collision in chain in- 
flation. As we will see, chain inflation, like any other 
inflationary model, additionally produces a red tensor 
spectrum at large scales through super-horizon quantum 
fluctuations of the metric. 

The most important contributions to the GW arise to- 
wards the very end of inflation, during the last few PTs. 
The peak frequency of the total spectrum is at /jv-i with 
amplitude £lhubh 2 (fN-i)- The contribution from much 
earlier PTs -which have a smaller peak frequency- are 
(nearly) invisible, because they are screened by the low- 
frequency tail of the spectrum of subsequent PTs due to 
the redshift factor e — 4 ( JV - m -- 1 )x, That is why the number 
of minima, N, does not change the final results as long 
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FIG. 3: The spectrum of GWs produced from chain inflation 
for x = 1/3 an d for various values of e. For lOGeV < e 1/4 < 
10 9 ' 5 GeV the signal is observable in some of the future GW 
probes. 



as N\ ^> I- However closer to the last bouts of chain 
inflation the effect of redshift factor diminishes; thus the 
effect of PTs that happened toward the end of chain in- 
flation are in fact observable in the spectrum as wiggles 

a t / ^5 /peak- 

Focusing on the frequencies around the tip of the spec- 
trum for chain inflation reveals a richer structure: wiggles 
in the GW signal due to different PT giving contributions 
at different frequencies. Since the resolution of Fig. 1 is 
inadequate to show these wiggles, Fig. © illustrates the 
tensor spectral index nx = dlnVlh^h 2 /2irdhi f close to 
the peak frequency /jv— i< This signature of wiggles near 
the peak from chain inflation can be used to distinguish it 
from an arbitrary single first order phase transition that 
happens to have the same peak frequency and amplitude 
at the peak frequency. Such wiggles can be a "smoking 
gun" for chain inflation. However, we caution that the 
details of the wiggles have not yet been computed pre- 
cisely, e.g. we have not considered the effects of having 
the m + 1th transition start before the mth is complete. 
Further work will be required to nail down the details. 

There is a second way to differentiate multiple vs. sin- 
gle tunneling events; i.e. to differentiate chain inflation 
from a single first order PT with the same peak frequency 
and amplitude. For both models to produce the same 
GW, the potential height Vq of a single PT would have 
to be related to the energy difference between vacua e of 
chain inflation in the following way: Vq ~ e/(l — exp _x ) 4 
(where x is the number of e-folds per PT in chain in- 
flation). The difference in the overall amplitude of the 
spectrum can be seen (with poor resolution) in Fig.l: 
matching both spectra at f max , one finds that the single 
PT has an overall amplitude that is lower than that of 
chain inflationl by up to 3% at frequencies below the max 
/ < /jv-i- The reason for this difference is the gradual 
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FIG. 4: Total GW spectrum produced during chain inflation 
from both bubble collisions and quantum fluctuations for \ = 
1/3, and for e 1/4 = 100 GeV and e 1/4 = 10 s GeV. 



decrease of the tensor spectral index at frequencies close 
to but smaller than /jv-i- This small difference in the 
amplitude at the low frequency end can in principle be 
detected. 

Three parameters in chain inflation determine the GW. 
The peak frequency in Eq (| 15[) depends on x (the number 
of e-folds per stage) , e (energy difference between any two 
vacua) , and (very weakly) TV (the total number of min- 
ima); while the GW energy density in Eq (fT5|) depends 
on x an d (again very weakly) N. Since the peak fre- 
quency and energy density scale exponentially with ~xN 
where xN >> 1, the iV-dependence can be ignored. The 
largest GW amplitude is found for the largest allowed 
value of Xi i- e -j X = Xc — 1/3. For the remainder of 
this paragraph we investigate this case and vary e. The 
GW amplitude remains the same for all values of e, and 
the spectrum just shifts to a different frequency. For 
10 7 - 5 GeV < e 1 / 4 < 10 9 5 GeV, at least a part of the 
spectrum from bubble collisions falls into the advanced 
LIGO sensitivity region. For values of 10 5 2 GeV < e < 
10 7 5 GeV, the peak frequency shifts to the BBO region. 
Finally, for 10 15 GeV < e 1 / 4 < 10 5 2 GeV, the peak 
frequency shifts to LISA region. The gravitational wave 
signature of chain inflation from bubble collision cannot 
be detected in these experiments for e 1 / 4 < 10 GeV and 
X = 1/3. Alternatively, if we follow the results of [ttI ]. 
scales as low as e 1 / 4 ~ few O.OlGeV for x = 1/3 should 
be detectable. For X — \ an d e 1 / 4 < few eV, the am- 
plitude of the spectrum from bubble collisions on large 
(~ horizon) scales becomes bigger than the observational 
bound from WMAP, Q G wh 2 < 1CT 15 . This puts a lower 
bound on the value of e 1 / 4 > 10 eV, assuming that the 
nucleation rate takes the critical value, 7 = 

One can investigate the dependence of the GW on the 
value of x- We assumed so far that the nucleation rate is 
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at the critical value so that x = 1/3- For smaller values 
of Xi the peak frequency shifts to higher values whereas 
the amplitude of gravitational waves decreases. Since 
the peak frequency of the spectrum occurs at /jv— 1, for a 
fixed value of e, the peak frequency varies like x ^^ e - x ^ ■ 
The dependence of the amplitude of the spectrum on 
7 is more involved, but as the amplitude at the peak 
frequency is few x r2b^ 2 (/jv-i) 7 one would expect that 



the amplitude would depend on % as x 2 (l 



One 



should note that for values of x that do not satisfy the 
condition Nx 3> 1, the dependence of the peak frequency 
and its amplitude on x are respectively given as 



X(l-e-*) 



and X 2 (l-e-*) 

Spectrum of GWs from Quantum Fluctuations: 

In addition to GW from bubble collisions, chain inflation, 
like any other inflationary model, can produce a tensor 
spectrum from quantum fluctuations of the space-time. 
The amplitude of such quantum fluctuations is propor- 
tional to the overall amplitude of the potential, V: 



P T {k) 
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where the expression should be calculated at the moment 
of horizon-crossing. The relevant scale for comparison 
with CMB experiments corresponds to ~ N e e-foldings 
before the end of inflation, where N e is the minimum 
number of e-foldings required to solve the horizon prob- 
lem (N e ~ 60 for GUT scale inflation); we take the cor- 
responding potential height to be V*. QqIi 2 can be re- 
lated to the amplitude of primordial tensor perturbations 
through 0: 



n Q h 2 



AgwPt = A GW rP s , 



where A GW = 2.74 x l(r 6 (^) 1/3 and P s = 2.45x10 

at k„ = 0.002Mpc~ 1 and WMAP+BAO+SN constrain 
r < 0.22 (95%)C.L. [23]. In chain inflation, the slope of 
the tensor power spectrum, tit relates to r through the 
relation 
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which is different from the well-known consistency rela- 
tion r = —8ht for standard single-field slow-roll inflation 
(n.b., a different setup that violates the consistency re- 
lation is described in [lij]). The extra factor of 3 in the 
denominator of the right hand side (RHS) is due to the 
fact that in chain inflation, the scalar spectrum is [Io| 
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where s' is the equivalent of the slow-roll parameter as 
defined in 10]. Thus for a given r, the tensor spectrum 
of chain inflation which arises from quantum fluctuations 
is redder with respect to its slow- roll counterpart. For 



example, for r = 0.01, one finds Oq/i 2 = 6.58 x 10 17 
at = 0.002MPC- 1 and n T = -3.75 x 10~ 3 for chain 
inflation. The total gravity spectrum from chain inflation 
will be 



Qtoth 2 = ilbubh 2 + Qq}i 



(26) 



For values of r detectable in the coming decades, the 
first term on the RHS dominates on large scales while 
the second term dominates on small scales, as can be 
seen in fig.(0]). The QUIET experiment [25[ or the pro- 
posed CMBPOL [7J have the capability of detecting ten- 
sor modes down to at least r = 0.01 on Hubble scales, 
k = 0.002 Mpc" 1 (/ - 3.09 x 10~ 18 Hz). 

The scale of the transition from dominance of 
quantum-induced GW to bubble-induced GW depends 
on the values of K (or, equivalently, r), e and x- F° r 
(r = 0.01, e = 10 8 GeV, and x = §)> the transition oc- 
curs at k t ~ 2.5 x 10 12 Mpc" 1 (/ ~ 10~ 2 ' 4 Hz). Alterna- 
tively, for (r = 0.01, e = 10 2 GeV, and x = 3), the tran- 
sition takes place at k t ~ 4 x 10 6 Mpc^ 1 ; see fig. (U]). For 
e = 10 2 GeV, and x = \i quantum fluctuations will al- 
ways dominate on large scales as long as r > 1.39 x 10~ 30 
(i.e. V* > 8.73 x 10 8 GeV) [13 . For a given x, the 
peak frequency shifts to lower frequencies as e decreases, 
whereas the amplitude of the peak frequency remains al- 
most constant. For x = \ an d e 1 ' 4 < few eV, the am- 
plitude of the spectrum from bubble collision becomes 
bigger than the observational bound, f^cw^ 2 < 10~ 15 . 
This puts a lower bound on the value of e 1 / 4 > 10 eV, 
assuming that the nucleation rate takes the critical value, 

Conclusion: In this letter, we calculated the profile of 
GWs generated during chain inflation. Two mechanisms 
are responsible for generating the signal; bubble collisions 
(which dominate on small scales) and quantum fluctua- 
tions (which can dominate on large scales). We found 
three observational signatures for chain inflation. First, 
the signal could become detectable in upcoming large fre- 
quency probes of gravity waves: Advanced LIGO, LISA, 
or BBO. In fact, wiggles in the GW spectrum due to 
multiple phase transitions may be a "smoking gun" at 
frequencies below the peak. Second, the amplitude below 
the peak would be different from that of a single tunnel- 
ing model. Third, the consistency relation between scalar 
and tensor modes from quantum fluctuations is different 
from the standard one for slow-roll inflation and may be 
tested in CMB experiments. 

We also wish to comment on the connection with the 
landscape in string theory. If inflation happened through 
a series of rapid tunnelings in the string landscape, this 
work shows that future high frequency probes of GW can 
shed light on the structure of the landscape. Another 
possibility is that the many vacua in string theory may 
be connected by a combination of tunneling and rolling, 
alternating with one another. This case may generate the 
same GW discussed in this paper due to bubble collisions, 
as long as several of the last e-folds of inflation are due 
to tunneling. If there is rolling near the beginning of 
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inflation but tunneling near the end, then the quantum 
fluctuations would produce the standard result while the 
bubble collisions would be as described in this paper. A 
full analysis of the possibilities of alternating rolling and 
tunneling warrants another paper. 
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